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SECTION  1 


INTRODUCTION 

A problem  of  continuing  and  widespread  interest  In  electro- 
magnetic theory  is  that  of  wire  configurations  under  certain  speci- 
fic conditions.  The  transient  and  harmonic  analyses  of  thin-wire 
structures  have  been  the  subject  of  investigations  for  a number  of 
years.  An  early  treatment  of  thin  wires  was  given  by  Oseen  [1]  who 
applied  the  method  of  retarded  potentials  to  straight  thin  wires  and 
calculated  the  induced  current  on  a wire  by  a transient  incident 
wave.  Hallen  [2]  used  a slightly  different  form  and  derived  a pure 
integral  equation  for  the  induced  current  on  a thin  wire  and  used  it 
to  derive  analytical  expression  for  the  natural  frequencies  and 
current  distributions  of  the  natural  modes.  Since  then  many  investi- 
gations have  studied  various  thin-wire  problems. 

In  recent  times  Tesche  [3]  analyzed  the  thin  wire  scatterer 
! rom  the  singularity  expansion  point  of  view.  Wilton  and  Umashanker 
[s]  conducted  a parametric  study  of  an  L shaped  wire  using  the  sing- 
ularity expansion  method  (SFM).  The  EMI1*  interaction  with  a thin 
wire  above  a ground  plane  using  SEM  was  investigated  numerical Iv  by 
Shumpert  [5],  However,  in  more  recent  times  Crow,  et  al.  [b]  have 
conducted  an  in  depth  study  of  the  crossed  wire  structure  using  SEM. 

Even  though  many  of  these  papers  are  helpful  in  the  determina- 
tion of  induced  currents  and  charge  densities  on  thin  wire  config- 
urations, most  of  them  utilize  very  complex  and  elaborate  numerical 
* 

The  electromagnetic  pulse  generated  by  a nuclear  denotation  is 
generally  referred  to  as  the  EMP. 
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techniques.  Being  motivated  by  some  techniques  used  earlier,  the 
present  investigation  utilizes  primarily  analytical  techniques  for 
the  treatment  of  wire  configurations  in  the  proximity  of  a lossy 
ground . 

The  second  chapter  of  this  report  presents  the  time  harmonic 
analysis  and  considers  the  plane  wave  illumination  of  a single  wire 
over  a perfect  ground.  The  approach  used  follows  Taylor,  et  al.  ['] 
for  the  current  distribution  induced  on  a transmission  line  by  a 
nonuniform  field.  Subsequently  the  plane  wave  illumination  of 
crossed  wire  configuration  over  perfect  ground  is  considered.  The 
single  wire  and  crossed  wire  conf igurat ions  are  then  considered  to 
be  in  the  proximity  of  a lossy  ground  plane.  The  third  chapter  in- 
cludes transient  analysis  with  brief  description  of  SEN.  It  presents 
the  technique  for  determining  the  natural  modes  and  frequencies  for 
perfect  ground,  as  well  as  lossy  ground  considerations  utilizing  SE>1. 
Numerical  results  obtained  for  botli  single  wire  and  crossed  wire  con- 
figuration are  presented  in  Chapter  IV.  Chapter  V contains  the  con- 
clusions and  comments. 

The  accurate  determination  of  the  natural  frequencies  and 
natural  current  modes  on  a mathematically  tractable  conf igurat ion 
(e.g.,  the  crossed  wire  conf igurat ion)  is  an  early  step  in  the  study 
of  electromagnet ic  field  interaction  with  an  aircraft  that  will  pro- 
ceed to  more  complex  conf igurations.  This  and  other  topics  are 
considered  in  some  detail  in  this  report. 
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SI  T I"  ION  1! 


T1MK  HARMON  1 0 ANALYSIS 


J.l  Plane  Wave  1 1 1 urn  i ua  tJL  on  of  Che  Single  Wire 

Over  Per leoC  Crouud 

Some  time  .i<o , Bates  and  Hawley  (Sj  presented  a simple  first- 
order  treatment  of  tin'  scattering  of  plane  waves  from  an  infinitely 
thin  wire  above  a ground  plane.  In  that  paper  the  effects  due  to  a 
loss v ground  plane  are  neglected  but  thev  were  later  considered  by 
Leonard  Schlessinger  l1-1].  In  more  recent  times  Crow,  et  al.  [r>|  and 
[10 1,  have  also  conducted  similar  studies  for  thin  wires  of  finite 
length  above  a perfect  ground  plane  applvlng  purelv  numerical  tech- 
n iques . 

The  present  invest igat ion  considers  (utilizing  simple  analyti- 
cal techniques!  plane  wave  illumination  of  a single  wire  above  a 
ground  plane.  Initially,  the  perfect  ground  case  is  presented, 
later  the  lossv  ground  plane  case  is  discussed.  In  both  cases 
transmission  line  theory  is  used,  which  requires  the  wire  height  to 
be  much  less  than  the  wave  length  and  the  wire  length. 

In  this  section  we  obtain  the  expression  for  a single  wire  over 
perfect  ground  illuminated  by  a plane  wave.  The  general  formulation 
for  the  current  distribution  along  the  thin  cylinder  (or  wire)  due 
to  normal  plane  wave  incidence  has  been  obtained  by  Taylor,  et  al. 
[7]  for  a two  wire  transmission  line  with  arbitrarv  termination 
impedances.  For  the  problem  considered  here  the  single  wire  with 
its  electromagnet ic  Images  forms  a two  wire  transmission  line  with 


Terminated  two-wire  transmission  line  excited 
by  a plane  wave  with  normal  incidence. 
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V is  the  shunt  admittance  per  unit  length  and  Z the  series  impe- 
dance per  unit  length. 


J . 1 Plane  Wave  Illmninat ion  e:  a Crosse, : h'i_tv  Con:  ijtur_.it  ion 

Over  a Perfect  Ground 

Several  investigators  have  considered  the  plane  wave  illumi- 
nation of  wire  configurations  over  perfect  ground  and  in  free  space. 
Shumpert  [“>]  has  also  analyzed  thick  wire  configurations.  Crow.  et. 
al.  [11],  formulated  the  problem  using  Pocklington  type  integro- 
differentiul  equations  and  studied  it  by  means  of  the  SE>t.  Numeri- 
cal techniques  were  used  to  determine  the  natural  frequencies  and 
associated  current  modes.  Analysis  of  crossed  wires  in  a plane  wave 
field  was  done  by  King  (12];  however,  recently  Crow,  et  al.  [5] 
utilized  SE>1  techniques  for  perpendicular  crossed  wires  over  a per- 
fectly conducting  ground  plane.  In  order  to  avoid  the  complicated 
numerical  techniques,  an  attempt  was  made  to  obtain  the  natural 
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frequencies  and  associated  current  modes  utilizing  transmission 
line  theory  in  conjunction  with  electromagnetic  superposition. 

Consider  the  wire  configuration  oriented  as  shown  in  Figure 
2.3.  Initially  the  current  distributions  are  obtained  considering 
excitation  along  the  S.  , .1,  or  £ elements  separately.  Then  the 
net  current  distribution  is  obtained  by  the  summation  of  the  fore- 
going results  via  superposition. 


Case  A.  Consider  the  excitation  along  the  i element  onlv. 
° a 

With  reference  to  Figures  2.3  and  2.4  the  current  on  the  C element 

a 

is  obtained  by  treating  the  configuration  shown  in  Figure  2.3. 

There  are  basically  two  junction  conditions  that  apply  to  this 
configuration.  (1)  The  Kirchhoff's  current  law  must  be  satisfied 
at  the  crossed  wire  junction  and  (2)  junction  voltages  of  the  wires 
are  all  equal,  hence  7.<  is  infinite  corresponding  to  the  open  end 
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and  where  Z and  Z e are  the  impedances  seen  looking  into  the 

U>  I 

transmission  line  formed  by  the  l and  i,  element  with  their 
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images,  from  the  junction,  Z j is  the  equivalent  impedance  formed 
considering  Z and  Z impedances  in  parallel.  Accordingly 
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Using  (2.13)  and  (2.15)  v to  Ids 


Due  to  the  different  coordinate  axes  being  in  opposite  directions 

for  t and  elements,  there  is  sign  difference  between  (2.18) 

and  (2.19)  and  their  corresponding  equations,  (2.10)  and  (2.11). 

The  current  induced  on  i element  is 
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Case  £.  Once  the  current  distribution  over  different  elements 
is  obtained  due  to  excitation  along  4^  and  elements  the 

total  currents  can  be  easily  determined  utilizing  superposition 
theorem.  It  should  be  noted  that  the  current  distribution  due  to 

the  plane  wave  excitation  along  the  4 elements  can  be  inferred 

0) 

from  the  foregoing  by  using  symmetry.  To  obtain  the  total  currents 
for  the  entire  structure,  results  for  Case  A are  added  to  the 
corresponding  results  for  Case  B to  yield. 
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l (z.J-.)  - I ( 0 , j ui ) ~ 

w s i n k i 


(2.24.) 


where 


I to, J- ) 


r i k ' ^ 

j E sin  k h 1 -7.  tan  — r + Z tan 
11  _ * i - a 


__  

a t ui 


k>f  ] 


l - l_l 


a t 


(2.25) 


Similarly  for  the  element 


sin  k(i  - z' ) 

IfCz’.J.,)  - 1,(0.  ja.) j— 

1 1 sin  k » - 


F sin  k h 

k 2 ".»ia  k 4.  ,sln  ***'"  + sin  k k *' 

Ct  f 


t 


t 


(2.  2o) 


where 


a t u> 


,k<f) 

L J.7  . ,,,  I 

k''a)“] 

l 2 1 "./an  ( 

2 I 

— 1 

!Z  Z.J 

a t 

(2.27) 

Lastly  for  the  » element 


s in  k ( l - z ) 

I * I (O.jio)  

a .J1  sin  k i 

a 


4 E sin  k h 
o 


+ kZ  sin  k i [sln  k(-';a>  + sin  k 0 "sin  k zl 
ca  a a a 


where 


k 'f  ) , k 0 

.]  4 E sin  k h | (Z  +2  Z ) tan  , — — - , +Z  tani-~- 
I (O.Jw)  - 0 L 1 1 • 1 *’  1 2 J 

a k[(Z+Z,)Z  + 2 Z Z e ] 

a t u>  at 

(2.29) 

Note  that  the  Junction  currents  given  by  (2.2“>),  (2.27)  .,nd 
(2.29)  satisfy  Kirchhoff's  current  law  at  the  crossed  wire  junction. 


In 


. 3 Current  Source  Kxci  tuition  of  Crossed  Wires  Over  a 
Perfect  Ground  Plane 

To  obtain  all  the  natural  current  modes  and  frequencies  it  is 
convenient  to  consider  the  crossed  wire  structure  to  be  driven  from 
a current  source  as  shown  in  Figures  2.5  and  2.6.  Then  applying 
transmission  line  theory,  the  current  and  voltage  along  the  i 
element  are 

VMy.Jw)  - V,e~Jkv  + VJky  . 12.30) 

I(y,ju))  - [V1e"Jky  - V:e,ky]  . (2.31) 

i. 

CU) 

where  the  following  end  conditions  must  be  satisfied: 


1 ('  ,Jw)  - — [V,e"jkiw. -V.e_JkCu)]  - I 

U)  U)  Cs  l - O 

Cu) 


(2.32) 


and 


V (0, j.j) 

u) 

I (O.ju.) 


V V2 

cu  V , - V , 


-1 

\~  + v-  + 1 . (2.33) 

L L L 

|_a  f u)  J 


On  performing  mathematical  manipulation,  (2.32)  and  (2.33)  yield 


, Z (Z  -Z  ) 
, ^ Cu)  Cu)  T 
J 1 


Z + Z_  Sin  k 

U)  I 


(2.34) 


, (Z„  + Z_) 

l Cui  c..  1 


+j  i z +t; 


$ in  k 


(2.  35) 


Using  (2.34)  and  (2.35''  in  (2.31)  yields 


i (y.J‘») 


i;[Zr  (e 

O Cu> 


-Jkv 


. Jky 


) - 


T 


(e 


-Iky 


-ejk>’)] 


(2.36'' 
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There  tore 


I (y.J-i) 

SJJ 


Cui 


cos  k v + 


sin  k v 


(I-  + Z_)  sin  k 

.0  r 


I (y.Jw) 


Cu> 


cos  k v + 


sin  k v 


i cos  k v + I Z,  sin  k v‘  o 

CIO  iO  ' 1 u> 


where 


a f 


Zf  Z.+  Za 


"a  f 


U\38) 


(2.  34) 


(2. AO) 


and  is  the  parallel  combination  of  all  elements  except  the  one  that 

is  excited.  The  current  on  the  element  that  is  not  excited 

u) 

and  on  the  i and  C . elements  mav  be  obtained  bv  the  use  of 
at 

Kirchhoff's  current  law  at  the  junction. 


-I  (O.jui)  - 1^(0, Jui)  + I'(O.jvo)  + lf(0,ju>)  (2.41) 

and  the  Junction  voltage  conditions, 

V (O.juO  ■ - I (0,Jui)Z  - -I.  (0,  jio)Z  - -I'lO.juOZ  (2. 42) 

UJ  a 3 I I u)  w) 

Hence  the  current  on  the  ^ element  that  is  not  driven  is 
Z Z \ 

+ l’(0,ju))  - -I  (0,ju)  . (2.43) 


:he  currents  on  other  elements  of  the  wire  configuration  can  be 


obtained  similarlv. 

They  are* 

sin  k(,  - v') 

i;(y',ju) 

a 

I'(O.jw)  — f 

sin  k v 

(2.44) 

,jw) 

a 

I (0,  ).„)  sin  k ( V z) 

a ; ■ 

S 111  K 

a 

(2.453 

a 

lff0,jw)  sin  k(vf  " z’  ] 

(2.463 

sin  k . 

t 

where  from  i2.43)  and  (2 

. 14) 

l’(0,ju) 

- + 

, T cj  1 o 

Z (Z  -r  2^)  sin  k > 

(2.47) 

Correspondingly  for 

the 

1 a and  f elements. 

Z_  Z,  I 

Ta(0.j  ) 

« + 

i A C U)  o 

2 (2  , + 2T  ) sin  k «: 

(2.48) 

7 7 T 

*-T  4 L 

^-(O.jw) 

= + 

. 1 c.i)  O 

Z . (Z  + Z ) sin  k ; 

1 U)  I a’ 

(2.49) 

When  the  element  is  excited  bv  a current  source,  the 

current  on  i , i.  and  v elements  can  be  derived  in  a similar 

cl  l U) 

manner  and  are  presented  as  follows 


Note  that  the  v axis  coincides  with  the  y axis  except  that  it 
is  directed  in  the  opposite  direction 
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| 


4 


:\ 


sin  k( 


fl 


l-  (V,J  A ' sin  k " 


sin  k(C  - z) 

Ia<z,.M  ■ Ia(0.Ju)  ,ln  kJ, 

a 


where 


Z Z . I 
c f o 


I'(0,Ju>)  = +J 


ZJZ{  + Z.)  sin  k .•  f 


I (O.Ju)  - +j 


ZV  Z . . I 
let  o 


Z (Z,.  + Z.)  sin  k L 

u)  t * 


la(0,ju)) 


= +j 


z,  z . 1 

let  o 


Z (Z  -f  Z ) sin  k v , 
at’.  t 


(2.59) 


(2.60) 


(2.61) 


(2.62) 


(2.63) 


and 


has  been  defined  previously. 


2 . 4 Lossy  Ground  Plane  Considerations 

The  propagation  constant  encountered  in  the  earlier  sections  is 
a real  quantity  for  the  perfect  ground  case  only.  However,  under 
lossy  ground  conditions  the  propagation  constant  becomes  complex. 
Although  only  a few  investigators  have  considered  the  treatment  of 
wires  over  a lossy  ground  plane,  there  has  been  remarkable  progress. 
Schlessinger  [9]  derived  a self  contained  expression  for  the  propa- 
gation constant  following  Sunde  [13]  and  Wait  [14]. 

Sunde  [13]  derived  the  propagation  constant  for  insulated 
aerial  conductors,  accounting  the  finite  conductivity  of  the  earth 
which  gives  rise  to  an  increase  in  the  longitudinal  impedance  and  to 


resistance  losses  in  the  earth.  He  derived  this  expression  as  a 
special  case  of  wire  of  infinite  length  above  the  surface  of  the 
earth. 

According  to  Sunde  the  series  distributed  inductance  can  be 
approximated  with  satisfactory  accuracy  lor  engineering  purposes' 


-’ll  + > h) 


L - - a 


( 2 . o4 ) 


whe  re 


= > juiu  (c  + juic  ) (2.65) 

with  and  c the  ground  conductivity  and  permittivity,  re- 

spectively. The  shunt  distributed  capacitance  is  approximated  by 


C 


2h  , 

,ni  - ) 


(2.66) 


According  to  the  foregoing  the  product  of  the  distributed 
impedance  and  shunt  admittance  is 

(Jw)-e  u 


ZY 


> n (— ) 
a 


.n;1+vi+  . 2h  n 

I ai  - •- — p-  + . n I 

> h : a i 

I 1 S J 


(2,67) 


Since  the  propagation  constant  k = -j>YZ  , then 


-J  >0  I 1 + 


. 2h 


1 + 


k n 


h) 

s 


- S 


> h 
g 


! i 


1 2 . pS) 


where  y = j ; 

o J o 0 


All  the  foregoing  expressions  for  the  current  induced  on  the 
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wire  configurations  over  a ground  plane  became  valid  for  the  los 
ground  when  the  propagation  constant  of  (2.68)  is  used  in  the 
transmission  line  current  and  voltage  expressions. 


Sli'TlON  111 


I'KANS  1 TNT  ANA1.YS1S 


1.1  Singularity  I'xyims  lyn  Mo  t hod 

The  STM  technique  was  rooontlv  introduced  bv  Baum  [IS)  as  a 
technique  tor  efficiently  characterizing  a conducting  bodies' 
t espouse  to  either  transient  or  steady  state  elect romagnel ie  (Mum- 
mat  ion.  An  arbitrarily  shaped  conducting  object  Is  completely 
characterized  bv  its  complex  natural  f requonc les  together  with  its 
associated  natural  current  mode  distributions.  The  mot  hod  mav  be 
viewed  as  an  extension  ot  ordinsrx  circuit  th.'orv  methods  to  distrib- 
uted parameter  systems.  Many  ot  the  recent  KMT  interaction  invest i- 
Sations  and  studies  have  utilized  the  S KM  technique.  Because  ot  its 
advantages  the  STM  technique  Is  applied  tv'  the  wire  configurations 
that  are  considered  In  the  foregoing. 

Initially  a l x 1 vector  is  formed  whose  elements  are  the 
current  distributions  on  the  crossed  wires.  Me., 


j V-*’] 

Ur. si  - : 1 l?.,s1 

, I 

■ 1 ly.sl 


where  s “ | and  Is  the  complex  frequency,  the  tilde  t 1 denv'tes 
the  l.aplace  transform. 

Then  utilizing  singularity  expansion  technique,  the  solut  ion 
tor  the  wire  currents  mav  be  expressed  as 


(3.  n 


Ur,s)  - l' 


o “ . s - s a 

it«  1 a 


l (r) 


n is  the  coupling  coefficient  and  i (r)  Is  the  ath  natural  mode 
and  where  the  natural  frequencies  s are  those  nontrivial  complex 
frequencies  obtained  by  the  complex  values  ot  s that  yield 


I.  (z'.s)  - 0 

t 


13.  J) 


I" si  - 0 


(3.31 


l ( v , s ) - 0 


(3.4) 


The  vector  i if)  represents  the  normalized  current  distribution 

l\ 

for  the  1 1 h natural  frequency  or  the  ath  natural  current  mode.  It 
Is  defined  as 


s - s 

t (r)  - =-  I ( r , s ) 

a s *s  c . r. 

a I o 


(3.5) 


where  c , Is  a constant  adjusted  so  that  the  maximum  value  of  any 
element  of  l^(,r)  (s  °qua ' to  1. 

Accordingly  then  the  coupling  coefficient  n ^ is  del  l nod  such 

that 


lm 


S “ s 

n 1 (r)  - — — — — l(r,s) 

a a s *S  t 

a o 


( 3 . h) 


However,  it  should  be  noted  that  the  natural  t requeue les  (or  simple 
pole  singularities)  of  llr.s)  occur  In  complex  conjugate  pairs  slnct 


tlu>  current  must  bo  purely  real  in  the  time  domain.  If  only  the 
frequencies  of  the  third  quadrant  o!  the  complex  s-plane  are  con- 
sidered then  the  time  domain  solution  for  impulse  excitation  is 


I ( r , t ) 


Ke  | n ^ oS.<(  i (,r1  | 


' . Natural  Frequencies  and  Associated  Natural  I'urrent 

Modes  for  a Perfect  tlround 

Plane  wave  Itlumination  ot  the  single  wire  over  perlect  g round 
has  been  formulated  in  section  .'.1.  In  this  sect  ion  we  obtain  the 
natural  frequencies  for  this  configuration  and  compare  the  results 
with  those  obtained  by  other  Investigators  utilizing  diltoront 
approaches . 

From  equation  l»’.l) 


■i  F s In  k h 

” - — j — r — | sin  ktr  - v 1 + sin  k » - sin  k z)  0.8^ 

k ..  sin  k > 


in  which  all  parameters  have  been  defined  earlier.  Note  that  the 


current  T(s,|,i')  diverges  at  natural  t requeue  i os  for  which 


sin  k v 


t l.  >M 


or  when 


t 1.  10> 


n - 1,  2,  l,  ...  where  v is  the  wire  length. 


For  the  perfect  ground  case  and  free  space  conditions  above  the 


ground  the  propagation  constant  is  purely  real 


k = - j / YZ  =■  ui*  LC  (3.  LI) 

where  the  characteristic  impedance  of  the  transmission  line  formed 
by  the  single  wire  and  its  image  is 


(3.12) 


Thus 


U) 


(3.13) 


relates  the  natural  frequencies  that  are  observed  to  be  independent 
of  the  height  above  the  ground.  These  natural  frequencies  being 
purely  real  differ  from  the  complex  ones  found  by  Tesche  [3]  for 
the  wires  in  free  space,  (a  different  physical  configuration),  but 
are  similar  to  the  ones  determined  by  Umashanker,  Shumpert  and  Wil- 
ton [16]  for  wires  over  a perfect  ground  plane. 

The  crossed  wire  conf iguration  is  more  difficult  to  analyze  in 
that  there  are  two  independent  conditions  for  which  the  induced 
currents  diverge  and  one  requires  a numerical  search  for  the  natural 
frequencies.  Both  divergence  conditions  occur  for  the  current 
source  excitation  (section  2.3)  but  onlv  one  condition  arises  when 
plane  wave  excitation  is  considered  (section  2.2). 

In  the  first  attempt  to  obtain  the  frequencies  for  which  the 
currents  in  (2.42),  (2.30),  (2.31)  and  (2.52)  may  diverge,  it  is 


r 


noted  that  both  (2.42)  and  (2.50)  diverge  as  Z 0 since 

U) 


7 2 

T V 0 - 


Z I 


(3.14) 


where 


i'(0,ju>)  — --  j -C1f  ° 
tu  7.  -*-0  2 sin  k Ji  Z 


cos  k l 

Z = -j  Z ■ 

w J cw  sin  k i 


Accordingly  the  natural  frequencies  are  the  frequencies  for  which 
Z = 0 . Hence  thev  satisfy 

U) 


ki  = (2a  - 1)  — = 1,  2,  3,  ... 


Or  in  the  complex  s-plane  the  natural  frequencies  are 


s = j (2a  -1)  — a = 1,  2,  3,  ... 


However,  it  should  be  noted  that  I and  1^.  remain  finite  at  the 

a f 

foregoing  natural  frequencies. 

A second  condition  for  which  the  crossed  wire  currents  diverge 


Z + Z_  =0 
co  1 


The  frequencies  for  which  the  foregoing  is  satisfied  can  be  obtained 
by  a numerical  search  routine  programmed  for  the  digital  computer. 


The  natural  current  modes  are  then  determined  for  the  associ- 


ated natural  frequencies  as  follows;  e.g.,  considering  the  first 


natural  frequency  kj> 


2.568. 

Bv  trial  and  error 

it  is  found 

and 

I at  k • k , . 

For  this 

i 

u>  1 

frequency  1 ^ is  then  adjusted  so  that 


l (Z.Ju)  ) 

a 1 


sin  k . ( >1  - z) 

* i.1 

sin  k,  i 
1 a 


If  (z'.Juij)  - -0.1506 


sin  k , ( - z 1 ) 
sin  k ( t 


I (y,ju)) 

J’l 


-1  (v'.uO 
10! 


-0.4291 


sin  k,  (C  - v' ) 
1 10 

sin  k , i 


(3.15) 


(3.1b) 


(3.19) 


Plots  of  the  first  six  natural  modes  for  crossed  wires  over  a ground 
plane  are  shown  in  Figures  4.1  through  Figure  4.6.  Analytic  results 
are  also  presented  in  Appendix  A. 

3. 3 Natural  Frequencies  for  a_  Lossy  Ground  Plane 

The  natural  frequencies  for  wires  over  a lossy  ground  plane  are 
obtained  by  using  Chapter  II  and  the  divergence  condition  (3.10). 
Hence  the  natural  frequencies  for  a single  wire  are  obtained  by 
finding  the  frequencies  for  which  the  following  is  satisfied 


1 4. 

1 

. in 

rvr 

| 2 

an 

1 ' 

• f2hl 

1 h 

C 

'iP  — 

a J 

l 8 JJ 

This  expression  must  be  solved  numerically  to  obtain  the  natural 
frequency  for  each  integer  a . Results  are  presented  in  Chapter  IV. 
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SUCTION  IV 


Nl'MF.RICAl.  RKSVT.TS 

The  numerical  technique  used  to  obtain  the  natural  frequencies 
for  single  wire  and  crossed  wires,  essentially  utilizes  search  rou- 
tines which  scan  the  area  of  interest  in  the  complex  s-plane.  Thus, 
the  natural  frequencies  which  satisfy  the  expression  in  question  are 
determined.  The  computer  program  utilized  for  t he  same  is  refer- 
enced by  Crow,  et  al.  [r].  It  uses  a subroutine  which  is  based  on 
the  Cauchy  Integral  formula  and  determines  the  zeros  of  the  analyt- 
ical function  I'izl  for  the  induced  current  expressions.  On  ob- 
taining the  results,  R expresses  the  ratio  of  the  function  F(rl 
to  the  average  magnitude  of  Fi.z)  evaluated  at  points  along  a se- 
lected search  contour.  It  should  be  noted  that  an  accurate  zero 
of  the  function  would  be  represented  by  values  of  R < 0.1 

Tables  4.1  and  4 . 2 show  the  natural  frequencies  for  a thin  wire 
above  a ground  plane  with  typical  conduct ivi t ies  and  dielectric  con- 
stants for  concrete  with  a » 0 . S m,  h « 1.0  in  and  C « 10.0  m 
(for  Table  -4.1)  and  a ” 2.0  m,  h ■ -j.O  m and  C » 40.0  m for 
Table  •< . 2 respectively.  In  Tables  ■> . and  -».-i  are  presented  natural 
frequencies  for  t ho  thin  wire  above  a ground  plane  for  different 

heights  with  conductivity  o ■ 0.0076  and  r ' ■ dielectric  con- 

81  8 

stant  ■ 25.0  for  a - 0.5  m and  - 10.0  m.  Table  *.0  is  valid 
for  n ■ 1.0  and  Table  4.4  is  valid  for  the  n “ 2.0  case.  The 
natural  frequencies  for  the  crossed  wire  st  rue ture  over  a perfect  ground 
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TABLE  4.3*  NATURAL  FREQUENCIES  FOR  THIN  WIRE  AT 
DIFFERENT  HEIGHTS  ABOVE  LOSSY  GROUND 
PLANE  (a  - 1 model 


a * radius  of  the 

wire  - 0.5 

m . v » length  of  the  wire 

“ 10.0m 

* ground 

c endue 

t i v i t v - 0 

,0076  mho/m. 

c * ■ dielectric 
a 

constant  « 2 

5.0 

U 

• c'  c ) 

S £ ° 

Height 

(h) 

(Sa/c) 

above  ground 
im) 

Rat  io 

(R) 

Real 

Imaginary 

1.0 

0.1001 

x 1 0-J> 

-0.05102-*14 

0.  2 7835b74 

1.5 

0.2331 

x 10" 5 

— 0 . 0307784*4 

0. 29960148 

2.0 

0.5121 

x 10-" 

-0 . 020^ 3dP4 

0. 30b099J0 

2.5 

0.b837 

x IO-" 

-0.01550837 

0 . 3089  J Jsb 

3.0 

0.1127 

x 10-  ’ 

-0.01215137 

0. 31 Os  18 38 

3.5 

0.  Ib5b 

x 10'" 

-0.00990076 

0.31131393 

4.0 

0.3717 

x 10-  ; 

-0.  DOS  U)  101 

0. 3119000b 

4.5 

0. 3b 2 7 

x IO"" 

-0.0071128 1 

0. 31230273 

5.0 

0. 2 1 bb 

x 10'" 

-0.0061997 

0. 31259685 

5.5 

0. 3973 

x ur" 

-0.0054785 

0. 31281843 

b.O 

0.4237 

x 10-" 

-0.00489o2 

0. 31299033 

b.  5 

0.8557 

x 10"' 

*•0  • 00m**  l 3 

0.31312b91 

7.0 

0.2168 

x 10'" 

-0.00401720 

0. 31323763 

7.5 

0. b56l 

x 10-’ 

-0.00 3b 7852 

0. 31332891 

8.0 

0. b5bl 

x 10'" 

-0.00338850 

0. 31 340527 

8.  5 

0.2344 

x 10-" 

-0.00313765- 

0 . .>  1 > *4  (■>  ^ 0 ^ 

9.0 

0.6883 

x 10'" 

-0.00391877 

0. 31352531 

9.5 

0. bbb9 

x 10'" 

-0.00272628 

0. 31357319 

10.0 

0.4135 

x 10- 

-0. 0025580b 

0. 3l3bls9s 

* With  reference  to  Figure  2.1. 


TABLE  4.4*  NATURAL  FREQUENCIES  FOR  THIN  WIRE  AT 
DIFFERENT  HEIGHTS  ABOVE  LOSSY  GROUND 
PLANE  (a  * mode) 


a - radius  of  the 

wire  » 0 

• 5m,  I 

» length  ot  the  wi 

re  » 10.0  m 

o * ground  conduc 

tivity  = 

0 . 00  7 o 

mho  m.  c 1 » dlelec 

trie 

constant  ■ 25.0 

U 

! 

K * 

* Js 

) 

Height  (h) 

(8 

/cl 

above  ground 
(ml 

Rat  iv 

.■>  (R'l 

Real 

I ma g i na r y 

1.2 

0. 1045 

\ 10" 

-0. 0*494  7 285 

0. ol*70751 

1.5 

0.  1423 

x 10"' 

-0.035*9272 

0. o201o536 

2.0 

0. 1050 

x 10" 

-0.0232108* 

0. 0238932 

2.5 

0. 1050 

x 10“- 

-0. 01t'47 lob 

0. o254  7 12 

3.0 

0.417-1 

x 10- 

-0.01307539 

0. o262959 

3.5 

0.  35 lt> 

x 10" 

-0.01057798 

0.  o2o78ol 

-..0 

0.2020 

x 10- 

-0. 00882*79 

0.6271034 

4 . 5 

0. 5290 

x 10- 

-0.00753400 

0.6273233 

5.0 

0.2724 

x 10-' 

-0.006548* 

0.6274822 

5.5 

0 • 1744 

x 10- 

-0.00577389 

0.6 2 760 17 

o .0 

0.1720 

x 1 0-v 

-0. 0031 5098 

0. p27o4*2 

o.  S 

0. 3420 

x 10-' 

-0. 00* c *030 

0. o277o7o 

7.0 

0.  1 74t> 

x 10- 

-0.00421481 

0. 6278270 

7.5 

0. 3338 

x 10— 

-0.00385542 

0. 0278759 

8.0 

0.1723 

x 10— 

-0.00354826 

0. 6274168 

3. 5 

0.5032 

x 10— 

-0.003233 01 

0.6274514 

4.0 

0.3131 

x 10- 

-0.00305189 

0.6279810 

o . 5 

0. 1887 

x 10— 

-0.00284884 

0. o 280066 

10.0 

0.17*2 

x 10“' 

-0.00266932 

0.6280289 

*W  1 1 h reference  to  Figure  2.1 
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plane  are  shown  In  Table  ->.5  for  + >'  . » 1.0  L,  - ./l  - 0.5 

at  ta 

and  a ■ 0.05 L (where  L is  a scale  factor).  The  natural  current 
modes  for  these  frequencies  ttirst  through  the  sixth)  are  shown  in 
Figure  -..l  through  Figure  •..<•>.  A comparison  of  these  nodes  with  t h. 
numerically  obtained  modes  of  Crow,  et  al.  [b]  exhibits  very  close 
agreement.  Lastlv,  the  natural  frequencies  as  obtained  oy  Crow, 
et  al.  [ r> ) for  the  crossed  wire  conf igurat ion  accounting  the  radia- 
tion losses  (with  the  same  dimensions  as  considered  for  results  in 
Table  ■*.!')  is  shown  in  Table  s.b  for  the  purpose  of  comparison.  It 
should  be  pointed  out  that  radiation  losses  could  be  incorporated 
into  the  presented  analysis  using  the  procedure  suggested  by  Marin 
[17],  This  is  reserved  for  future  study. 

Coupling  coefficient  calculations  are  presented  in  Appenix  A 
which  completes  the  SIM  analysis  for  the  crossed  wires  over  a lossv 
ground  plane. 


TABLE  4 . 5*  NATURAL  FREQUENCIES  FOR  CROSSED  WIRE 
CONFIGURATION  OVER  A GROUND  PLANE  FOR 


;.  + ■.  = l , ; / ( , » l . 

at  a f 

source  excitation). 

9-*12  (using  current 

a 

S /Lc 
a 

i 

j 

2. 5680 

■> 

j 

3. 1416 

3 

j 

4 . 0157 

4 

j 

6.2S34 

5 

. j 

7.9421 

o 

j 

9.424S 

7 

j 10.9614 

8 

jl 

2.5665 

9 

jl 

4 64o8 

10 

j 15. 70S0 

A 

With  reference  to  Figures  2.3  and  2. 

5 

TABLE  4.6 

NATURAL  FREQUENCIES  FOR  CROSSED  WIRES  (as 

obtained  by  Crow,  et  al.) 

FOR  c + ( * L , 

1 c . = 1 , h / L * 0 • 1 , 
a r 

:/L  = 0.05. 

a 

S L/c 
a 

1 

-0.0202 

j 2.4525 

•» 

-0.0426 

j 2. 7604 

3 

-0.0470 

j 3.9166 

4 

-0. 1775 

j 5.93S9 

5 

-0.3015 

j 7.7876 

6 

-0.32S9 

j 8.3861 

7 

-0. 5359 

j 10. 7325 

8 

—0 . 8424 

j 12.0186 

9 

-0 . 7667 

j 14 . 1852 

•> 


igure  4.2:  Natural  current  mode  for  a crossed  wire 
structure  over  a ground  plane,  trans- 
mission line  formulation  results  for 
^ = = la  + c-a/«.f  = l.Q412  (arrows 

indicate  directions  assumed  for  positive 
current ) . 
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i^ure  -t.o:  Natural  current  mode  for  a crossed  v 
structure  over  a ground  plane,  trans 
mission  line  formulation  results  for 
L " :\j  " "a  + • '' a ' 1 f = 1-9412  ( 

indicate  directions  assumed  for  nos i 
current  1 . 


SECTION  V 


CONCLUSION  AND  COMMENTS 

The  study  of  wire  configurations  in  the  proximity  of  lossy 
ground  is  useful  in  the  EMP  simulation  studies  of  aircraft  [18]. 

As  compared  to  other  more  elaborate,  complex  and  purely  numerical 
techniques  utilized  for  this  interaction  problem,  the  theoretical- 
numerical  approach  presented  here  is  much  simpler.  Besides,  the 
step-by-step  approach  enables  the  reader  to  obtain  a clear  under- 
standing of  the  interaction  problem. 

In  case  of  a single  wire  over  a perfectly  conducting  ground 
plane,  the  natural  frequencies,  as  shown  in  section  3.2,  are  inde- 
pendent of  the  height  above  the  ground.  It  should  also  be  noted 
that  in  this  case  the  imaginarv  parts  of  the  natural  frequencies 
are  close  to  those  of  the  wire  in  free  space.  For  the  single  wire 
above  lossv  ground  plane  the  natural  frequencies  for  different 
ground  parameters  are  tabulated  in  Chapter  IV.  Tables  4.1  and  4.2 
reveal  that  for  a tvpical  set  of  parameters  such  as  radius  of  wire, 
height  of  the  wire  above  a ground  plane,  and  length  of  the  wire, 
with  the  change  in  ground  conductivitv  and  dielectric  constant, 
there  is  not  much  change  in  the  resonating  natural  frequencies 
(zeros).  But  Tables  4.3  and  4.4  indicate  that  the  real  part  of  the 
zeros  becomes  more  negative  with  the  decrement  of  the  height  of  the 
single  wire  above  a ground  plane.  The  agreement  of  this  observa- 
t ton  with  the  results  obtained  by  other  investigators  verifies  the 
presented  analysis  of  thin  wire  configurations. 


*3 


It  should  also  be  noted  that  the  crossed  wire  configuration 
illuminated  by  a plane  wave  does  not  excite  all  the  natural  modes. 
However,  the  crossed  wire  configuration  excited  from  either  end  by 
a current  source  excitation  does  excite  all  the  natural  modes.  The 
determined  natural  frequencies  and  the  associated  natural  current 
modes  thus  obtained  are  in  good  agreement  with  those  obtained  bv 
Crow,  et  al.  [6]  (as  shown  in  Tables  4.5  and  4.6,  and  Figures  4.1 
through  4.6).  . . 
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A1TKNDIX  A 


— .1  . -.-.Mlf 


sin  k 


v < 


1 i v . 1 


sin  k 


For  tin-  n - - mod#.  . " ''  • 


l U.1-  ' ‘ 

i •* 


l .0 


sin  k i 
sin  k 


» li»- t > 


sin  k i ' 


_ ) 


l l'H> 


-0. 


sink 


sink  v > - v ' ' 


I lv\  .1-  > ■ -0.001 


s v n iv 


For  the  n * s no  Jo . - *■  ' 


sin  k, i' ^ 

X.U.J«S>  ’ 1,0  ~ sTn  k . ; " 


who  t o 


sin  k i i ( ~ * 3 

If(*\Jw5>  “ +0*27S  sin  k tf 


sin  k.  (v  _ v ' 

I (y\J*<0  * -°*t’35  — sin  \l~T~ 

Ul  k1  u) 


Lastly  for  the  n - f mode,  ^ 


o.AJsS  c/L 


IaU.J-,;>  ’ 0 

" 0 


sin  k (i 


who  r 


For  n 


mod*? 


A 

'k.  1 

•k;  J 

-0.001 

— 

= - ini 

j (k  - k 3 

--  f tan  + 

i * z f 

2 tan 

A 

t j 

k -k  , 

iiz.<z.,  - zf>  * 

J - 

a "f 

i 

k 

at  k = k 

+ 

Too  * 

where  k 

■ 6*284  I . , one 

obta ins 

c 

i r 

* 

0.0048 

n = 5 mode 

* 

!-) 

Jli 

• k ) 
a 

2 tan  | 
a 

k . 

• 

-0.635 

s v im 

j (k  - k ) 

~ r wan  ^ j + 

•> 

I c J 


k[z  (z  + z.)  + 2 z z 

--it  a f 


Then  at  k - k + 10Q  , where  k?  = 7.9421/L  , one  obtains 


' c ^ 0.00079 


For  n = b mode 


i.o 


• un 

l.  c k-k. 


k k 

-Z  tan  l — • + Z tan  -~ 
j (k  - k , ) — 1 c- = a 

k[Z  (Z  + Z ) + :z  Z.] 
a t at 


Then  at  k = k + — — 

- 100 


. where  k * 9.42-*8/L  . one  obtains 
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